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Abstract. In this paper we study the ground state phase diagram of a one-dimensional t—J—U model away 
from half-filling. In the large-bandwidth limit and for ferromagnetic exchange with easy-plane anisotropy 
a phase with gapless charge and massive spin excitations, characterized by the coexistence of triplet super- 
conducting and spin density wave instabilities is realized in the ground state. With increasing ferromagnetic 
exchange transitions into a ferrometallic and then a spin gapped triplet superconducting phase take place. 

PACS. 71.10.Hf - 71.10.Fd - 74.20.Mn - 71.27.+a - 75.10.Pq 

1 Introduction [22]. This triggered interest in low-dimensional models of 

correlated electrons showing mechanism for Cooper pair- 
Experimental and theoretical investigations show that many ing of triplet symmetry coexisting or closely competing 



strongly correlated electronic systems exhibit surprisingly 
complex phase diagrams £Q. Especially the competition 
or coexistence of magnetic order and superconductivity is 
a topic of increased current interest. Magnetically medi- 
ated Cooper pairing near the antiferromagnetic instabil- 
ity is widely discussed in the context of superconductivity 
in copper-oxide systems Moreover, the discovery of 

Triplet Superconductivity (TS) in S^RuO^ 0] and the 
recent discovery of coexistence of the TS phase with fer- 
romagnetism in U Ge 2 , U RhGe jSj and ZrZn 2 [7| has 
triggered an increased activity in studies of correlated elec- 
tron models showing close proximity of triplet supercon- 
ducting and ferromagnctically ordered phases [STITJlllUlllH 
I12JIi3U14II15II161I17| . 

Another challenging problem is superconductivity in 
quasi-one-dimensional compounds |19| . More than two decades 
have passed since the discovery of superconductivity in so- 
called Bechgaard salts such as (TMTSF) 2 X with X=PF 6 , 
CIO4, etc. [20]. Most interesting is the phase diagram of 
(TMTSF) 2 PF 6 which shows a spin-Pcierls (SP) phase 
in the ground state at atmospheric pressure. Increasing 
pressure leads first to a transition from the SP phase into 
a spin density wave (SDW) phase, and finally to the sup- 
pression of the SDW ground state in favor of superconduc- 
tivity |21| . In last years an increasing amount of convincing 
experimental evidence has been accumulated in favour of 
the triplet superconducting ordering in these compounds 



with SDW type ordering. Although a wide variety of the- 
oretical approaches have so far been used to study this 
topic ^23 , 24 , 25 , 26 , 27 , 28 , 29] , the understanding of the mi- 



croscopic mechanisms for the complex phenomena in these 
compounds is still far from complete. 

In an extension of the Hubbard model includ- 
ing anisotropic spin-spin interactions has been proposed 
as a suitable model for systems with coexisting orders. 
The model describes a system of itinerant electrons with 
anisotropic spin-exchange interaction between electrons 
on nearest-neighbor sites. The one-dimensional version of 
the Hamiltonian reads: 
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Here c* n a (c„ a ) is the creation (annihilation) operator 
for an electron at site n with spin a, p. a (n) = c? n a c n a , 

S(n) — \(^ n a a ' a pC n a where a 1 (i — x,y,z) are the Pauli 
matrices. 

Indeed, this model was shown to exhibit an extremely 
rich phase diagram in the case of a half-filled band PHI 
131) . In particular in the case of the ferromagnetic XY- 
typc exchange interactions the ground state phase dia- 
gram consists of sequence of transitions (with increasing 
ferro exchange) from a metallic phase with coexisting TS° 
and SDW Z instabilities, into an insulating Neel type anti- 
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ferromagnetic phase and finally, for strong XY-type ferro- 
exchange into the insulating ferromagnetic XY phase |31| . 

In this paper we study the effect of doping on the 
ground state phase diagram of the model 0). We use a 
combined approach based on continuum-limit bosoniza- 
tion and DMRG techniques. In the range of applicabil- 
ity of the continuum-limit approach we have obtained the 
ground state phase diagram for a wide range of model pa- 
rameters and band-fillings v. In our numerical studies we 
restrict our consideration to the case of a quarter-filled 
band with XY-type spin exchange interaction and on-site 
repulsion U > 0. We investigate the excitation spectrum 
of the system as well as the behavior of various correlation 
functions. Depending on the relation between the model 
parameters J xy /t and U jt we have shown the existence 
of four different metallic phases in the ground state (see 
Fig. II Op t : In the case of antifcrromagnetic exchange the 
system shows properties of a spin gapped (Luther-Emery) 
metal with coexistence of the singlet-superconducting and 
charge-density-wave (CDW) instabilities. The line J xy = 
corresponds to the Luttinger Liquid phase and marks the 
transition from a spin gapped sector with singlet Cooper 
pairing (antifcrromagnetic exchange) into the spin gapped 
sector with triplet Cooper pairing (ferromagnetic exchange) . 
For weak ferromagnetic exchange, at J° y < J xy < the 
system displays properties of a spin gapped metal with co- 
existence of the triplet superconducting and spin-density- 
wave (SDW Z ) instabilities. For ferromagnetic spin cou- 
pling the spin gap dependence on strength of the trans- 
verse exchange exhibits a dome-type shape, opening at 
Jxy = and closing at J xy = J xy . At J xy < J xy a rather 
unconventional fcrromctallic phase with gapless spin exci- 
tations and strongly dominating transverse ferromagnetic 
and triplet superconducting instabilities in the ground 
state. Finally, at J xy < J}~, our numerical data indicates 
on opening of the spin gap. In this phase the system is 
expected to show properties of a triplet superconductor. 

The paper is organized as follows: in the next section 
the weak-coupling continuum-limit version of the model 
is investigated. This allows to derive the weak-coupling 
phase diagram (Sec.|3J). In Sec. ^results of DMRG studies 
for chains up to L — 120 sites are presented. Finally, Sec. [3] 
is devoted to a discussion and concluding remarks. 



2 Continuum-Limit Bosonization 

Below in this section we consider the low-energy effective 
field theory of the model in the case of non-half-filled 
band. 

The standard bosonization procedure allows to express 
the initial lattice model in terms of two independent bosonic 
Hamiltonians 

TL = 7i c + Ti. s 

describing respectively the charge (c) and spin (s) de- 
grees of freedom. For the band-filling v ^ 1/2, the gapless 
charge sector is described by the free Bose field Hamilto- 
nian 



J dx{-^{d x ip c ) 2 + K c (d x tf c ) 2 }, (2) 



while the spin sector is governed by the quantum Sine- 
Gordon field 

H, = y J dx[ -^(d x <p S f + K s {d x $ s ) 2 

i(V^ips)}- (3) 



2m, 



Here ip CiS (x) and $ c . s (x) arc mutually dual bosonic fields 



dtVc 



[(p c ^(x),-d c . s {x)] 



(4) 



and ao is the infrared cutoff of the theory. The model 
parameters are given by 
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U - J x , 



— cos(27r;/) 
2 v ; 



(7) 



and the velocities of charge and spin excitations w c ,(s) = 
vf / K c ,(s)i where vp = 2tao sin(7ri/). 

Since at v ^ 1/2 the charge sector is described by the 
free Gaussian field J5| the vacuum averages of exponen- 
tials of the charge fields show a power-law decay at large 
distances 



/ g i \/2Tr Lp c ( X ) g - i \f2TTip c ( X ' ) \ 



-1/K C 



(8) 
(9) 



and the only parameter controlling contribution of the 
gapless charge degrees of freedom to the infrared prop- 
erties of the system is the charge LL parameter K c . 

The infrared behavior of the Sine-Gordon Hamiltonian 
Tls is described by the corresponding pair of renormaliza- 
tion group (RG) equations for the effective coupling con- 
stants K s (I) and M s (l) 



dM s (L) 

dL 
dK s (L) 

dL 



-2{K S (L) - 1)M S (L) 



(10) 



where L = ln(ao), K S (L = 0) = 1 + hgos an d M S (L = 
0) = g±/2n. The pair of RG equations l|10(l describes 
the Kostcrlitz-Thoulcss transition [33] in the spin chan- 
nel. The flow lines lie on the hyperbola 



Mi 



9os 



(11) 



and depending on the relation between the bare coupling 
constants go s and g± exhibit two different regimes [84| : 
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Weak- coupling regime. For go s > \g±_\ we are in the 
weak-coupling regime with effective mass M. s — ► 0. The 
low energy (large distance) behavior of the corresponding 
gapless mode is described by a free scalar field. 

The vacuum averages of exponentials of the correspond- 
ing fields show a power-law decay at large distances, 



where 



/ e iV2x<p s (x) e -iV2n<p s (x')\ ^ \x — x'\~ K " , 



(12) 
(13) 



and the only parameter controlling the infrared behavior 
in the gapless regime is the fixed-point value of the effec- 
tive coupling constants K* = K s (l = oo) determined from 
Eq. ((TTJ). 

Strong coupling regime. For 2(K® — 1) < \M®\ the 
system scales to strong coupling: depending on the sign of 
the bare mass M s °, the renormalized mass A4 S is driven to 
±oo, signaling a crossover to one of two strong coupling 
regimes with a dynamical generation of a commensurabil- 
ity gap in the excitation spectrum. The flow of \M S \ to 
large values indicates that the AdsCOsV&K^s term in the 
sine-Gordon model dominates the long-distance properties 
of the system. Depending on the sign of the mass term, 
the field ip s gets ordered with expectation values [12] 



S„ S sm (x) + S st {x) (18) 

S z sm {n) == O fm (x) ~ -L^.(s), (19) 

S* m (n) = O x FM {x) 

~ sin(\/27r^ s ) cos (%/2?n9 s ) , (20) 
SLW = ^ M (x) 

~ sin(v^7T^ s ) sin U/^hid^ , (21) 



and 







(M s ° > 0) 
(M s ° < 0) 



(14) 



S; t (n) = SiJW r.(x) 

~ cos(v / 27riy9c + 2k F x) sin (^/2~KLp^j , (22) 
S^(n) = Osdw-(x) 

~ cos(v / 27riy9 c + 2k F x) sin ^\/27n? s ^ , (23) 
5 s 2 t (n) = O sw ,(x) 

~ cos(%/2tt^ c + 2fc F a;) cos (-\/27n9 s ) . (24) 

In addition we use the following set of superconducting 
order parameters: 

3a) the on-site singlet 



It easy to check that, using the initial values of the 
coupling constants given in we obtain the following 

condition for generation of a gap in the spin excitation 
spectrum 



U - J xl 



cos(27w) 



3b) the extended singlet 



(25) 



> 



J. 



U — J xy cos(27w) — J z + — cos(27w). (15) 



3c) and the triplet pairing 



2.1 Order parameters 

To clarify the symmetry properties of the ground states of 
the system in different sectors we consider the following 
set of order parameters corresponding to the smooth "sm" 
and staggered "sf parts of : 

1) the on-site density operator 

p(n) => p sm (x) + p s t(x) 

where 



V2 



1 



c f c f 



(26) 



TSO (x)+0^ TS o(x) 



,T ^ C n,l C n+l 



O 



(X). 



(27) 
(28) 



Here operators without and with the subscipt rj corre- 
spond, respectively, to the smooth and 2fci?-modulated 
(staggered) parts of the corresponding superconducting 
order parameters. 

The bosonized expressions for the smooth parts of the 
corresponding superconducting order parameters, up to 
the accuracy of irrelevant ^-dependent amplitudes and 
phase shifts are given by 



p sm {x) ~ y -d x ip c (x), (16) 

p st {x) = C dw{x) 

~ sin(v / 27Ti J 9 c — 2k F x) cos(v / 2~7ri J 9 s ), (17) 

2) the on-site spin-density 



°ss( x ) ~ °es( x ) ~ cos(V2ntps) e 



/2tt& c 



o 



O\, s o{x) ~ sin(V27rv ;, s )e* 

\os(V2^^ s )e l ^- a " 
s\n(y/2^d s )e 1 ^ B -' 



TS^y 



(x) 



(29) 
(30) 

(31) 
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Similarly the bosonized expression for the staggered com- 
ponents of the corresponding superconducting order pa- 
rameters are given by 



sin(v / 27r<f3 c + 2k fx) e l 



(32) 

Note that the smooth part in Eq. I|25[) corresponds to 
the usual BCS-type pairing while at half-filling the oscil- 
lating terms in (|25|) and (|26|l describe the weak-coupling 
analogs of the //-pairing superconductivity |36| . 



3 The Weak-Coupling Phase Diagram 

In this section we consider the ground state phase dia- 
gram of the model Q away from the commensurate value 
of the 1/2-filled band. Due to the invariance of the model 
parameters under the transformation v — > 1 — v we re- 
strict our consideration to the sector < v < 1/2. Below 
we consider the the weak-coupling ground state phase at 
quarter-filling in detail, while for v ^ 1/4 we present a 
qualitative description of the phase diagram. 



TS°+ SDVV' 

( FM Z +il - Super conductivity ) 



Spin gap 

<cp >=K/2 •••• A 



SS + ES + CDW 

z 

( FM + T| - Super conductivity ) 

Spin gap 

«p s >=0 



B 




Fig. 1. The ground state phase diagram of the 1/4-filled itin- 
erant t — J xy — J z model. The thick solid lines indicate borders 
between the weak-coupling limit phases. Along these lines and 
in the sector C the spin excitations are gapless. The charge ex- 
citation spectrum is gapless in all sectors of the phase diagram. 



The longitudinal ferromagnetic and the ^-superconducting 
correlations decay faster as 



3.1 The U — case. 

Let us start with the case U = where the basic equations 
read: 



Kr 



1 



1 



2ttvf 



J xy 



1 



J z I cos(2-kv) . 



(33) 



and 



2(K S - 1) ~ J z + (J xy - i J z ) cos(2^) , (34) 



27rm s 



l 

TTVp 



Jxy + \ Jz C0S(2-KV) 



(35) 



In the following we study the two cases v = 1/4 and 
v 7^ 1 /4 separately. 



3.1.1 v = 1/4 

At U = the charge sector is featureless with K c = 1 in- 
dependently from the values of the exchange couplings J z 

cllicl Jx y • 

Since the charge sector is featureless, the ground 
state phase diagram is completely determined only by the 
spin degrees of freedom. The condition for dynamical gen- 
eration of a spin gap reads | — J xy \ > — J z and thus the spin 
sector is gapless for J z < — \J xy \ and along the semi-axis 
Jz > 0. These conditions determine the following three 
sectors of the phase diagram (see Fig. H): 

The sector A (J xy < 0,J Z > —J xy ) corresponds to 
a spin gapped phase with dominating SDW Z and TS° 
ordering. These correlations decay as power-laws at large 
distances: 



(0 8 DW>(P)0sDW>(r)) 
(0 TS o(0)0 TS o(r)) 



cos(7rr/2) 



(36) 
(37) 



(ew(o)ew(r)) = 

(O^ss(O)O^ss(r)) 



(38) 

- (O v -Es(0)O^Es(r)) 



~ (0,_TSo(0)0,_rso(r)) ^ cos(^r/2) • r~ 2 . (39) 

Other correlations decay exponentially. 

The sector B (J xy > 0,J Z > ~J X y) corresponds to 
a spin gapped phase with dominating SS, ES and CDW 
ordering. The corresponding correlations show a power- 
law decay at large distances: 



(O S s(0)Oss(r)) = (O E s(0)O ES (r)} ~ r" 



(40) 



(O C Dw(0)O C Dw(r)) £Si cos(7rr/2) • r" 1 . (41) 

In this case also the longitudinal ferromagnetic and the rj- 
superconducting correlations decay faster (|38|) and Ij39(l . 
while all other correlations are exponentially suppressed. 

In the sector C (J 2 < — \Jxy\) the gapless Luttingcr 
Liquid (LL) phase is realized. All correlations decay as 
power-laws, however the transverse antifcrromagnctic and 
triplet superconducting instabilities arc the dominating 
instabilities in this sector. The corresponding correlations 
decay as 



{OsDW*y {Q)OsDW*y i r )) 
(O TS *v(0)O T s*v(r)) 



cos(7rr/2)r 1 
-1-1/ K* 



(42) 
(43) 



where the fixed point value of the spin LL parameter 

1 



K* = 1 



2ttvf 



Jz Jx y > 1 • 



(44) 



Since the spin stiffness parameter K* > 1, the longitudi- 
nal ferromagnetic and the ry-superconducting correlations 
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given by Eqs. I|37(l and l|38(l decay faster and are the sub- 
leading instabilities in this sector. The transverse ferro- 
magnetic correlations decay even faster as 

(0 FM ~.v(P)0 F M-Ar)) - r- K W m ' . (45) 

However, within the accuracy of the used first order RG 
approach their decay is the same as of the longitudinal 
ferromagnetic correlations l|38l) . 

Finally the CDW, SDW Z and the TS° correlations 
decay as 

(OcDw(O)OcDw(r)) ~ {0' SDW (p)0' SDW {r)) 

~ cos(7rr/2) ■ r" 1 '^ , (46) 
(O TS0 (0)O TS a(r)) -r- 1 -^, (47) 

and correspond to the weakest instabilities in this sector. 




Fig. 2. The ground state phase diagram of the 1/3-filled itiner- 
ant t — J xy — J z model. The thick solid lines denote lines in the 
parameter space where the spin gap opens and indicate borders 
between the weak-coupling limit phases. The thin solid line cor- 
responds to the line K c — 1 and marks the transition from the 
sector of phase diagram with dominating density instabilities 
(K c < 1) into the sector with dominating superconducting or- 
dering (K c > 1). The dashed lines mark (qualitatively) the 
crossover areas between the Luttinger Liquid phases with dif- 
ferent sets of subleading instabilities. The phases A and B, as 
well as the D2 and E2 phases are characterized by an identical 
excitation spectrum and interchange of leading and sublead- 
ing instabilities. The subleading instabilities are indicated in 
brackets. 



3.1.2 v ^ 1/4 

At v 7^ 1/4 the charge stiffness parameter K c ^ 1 and 
therefore the line J z = —2J xy divides the parameter space 
into two semiplanes: the part with dominating CDW or 
SDW instabilities at K c < 1 and the part with dominat- 
ing superconducting instabilities at K c > 1. However, the 
effect of charge sector essentially depends on the band- 
filling. 



At 1/4 < v < 1/2, K c - 1 ~ sign(J z + 2J xy ). There- 
fore the dynamical generation of a spin gap and subse- 
quent pinning of the spin field with vacuum expectation 
value (<p s ) = ir/2 results to metallic phase with dominat- 
ing Triplet Superconducting (TS°) instability at J xy < 0, 
Jxy < Jz < —2 J X y and dominating antiferromagnetic 
SDW* at -0.5 J z < J xy < -0.5cos(2ttv)J z (see Fig. 
sectors A and B respectively). 

In the sector C of the phase diagram, at —2J xy cos 2 (-7w) < 
J z < —2J xy /cos(2'Kv), the spin gapped metallic phase 
with dominating CDW phase is realized. 

In the Luttinger Liquid sectors of the phase diagram 
D and E, at J z < mm.{J xy , —2 cos 2 (irv) J xy } the line J z = 
—2J xy marks the transition from a LL phase with dom- 
inating TS xy instability at J z < min{ J xy , — 2 J xy } (see 
Fig. [21 sector E) to a LL phase with dominating trans- 
verse antiferromagnetic instabilities SDW xy at —2J xy < 
J z < —2 cos 2 (m/) J xy (see Fig- El sector D). One has to 
note, that for a given band- width v ^ 1/4 in the LL sec- 
tors D and E, there are additional "subdominant order 
crossover" lines (marked by dashed lines in Fig. [21 and [3J 
which separate two areas with two different subdominant 
order. In particular, in the subsector E\ the subdominant 
order is SS + ES + TS° and in the subsector E 2 it is 
SDW xy . Similarly, in the subsector D\ the subdominant 
order is CDW + SDW Z and in the subsector D 2 it is 
TS xy . 




Fig. 3. The ground state phase diagram of the 1/6-filled itiner- 
ant t — J xy — J z model. Solid lines indicate borders between the 
weak-coupling limit phases. The thick solid lines denote lines 
in the parameter space where the spin gap opens and indicate 
borders between the weak-coupling limit phases. The thin solid 
line J z = —2J xy marks the transition from the sector of phase 
diagram with dominating density instabilities (K c < 1) into 
the sector with dominated superconducting ordering (K c > 1). 
The dashed lines mark (qualitatively) the crossover areas be- 
tween the Luttinger Liquid phases with different sets of sub- 
leading instabilities. The phases A and B, as well as the D2 
and E2 phases are characterized by an identical excitation spec- 
trum and interchange of leading and subleading instabilities. 
The subleading instabilities are indicated in brackets. 
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On the other hand, at v < 1/4 (sec Fig.^J the A sec- 
tor corresponds to the spin gapped metallic state with 
dominating SDW Z instabilities, while triplet supercon- 
ductivity TS° is the dominating instability in the nar- 
row stripe of sector B. Similarly, in the sector C at J z > 
max{— 2 J xy / cos(27ri/), — 2 J XJ/ cos 2 {-kv)} a spin gapped phase 
with dominating tendencies towards singlet superconduct- 
ing ordering is realized. Similarly in the Luttingcr liquid 
phase dominant and subdominant instabilities in the sec- 
tors D and E change places, i.e. in sector D the dominant 
order is TS xy and the subdominant order in sector D\ is 
SS + ES + TS° and in sector D 2 it is SDW xy . In sector 
E the dominant order is SDW xy and the subdominant 
order in sector E\ is SDW Z + CDW and in sector E 2 it 
is TS xy . 



3.2 The J z = case 

At J z = the basic equations read 
1 



1 - 



2ttvf 



(U - J xy cos(27n>)) , 



and 



2(K S - 1) ~ — !— U - J xy cos(27w) 
ttvf 

1 



27rm s 



TTVp 



U — J- 



•'■y 



(48) 

(49) 
(50) 



Again we distinguish the cases v = 1/4 and v ^ 1/4 in 
the following. 



3.2.1 v = 1/4 

At U 7^ 0, even for i/ = 1/4 the charge stiffness parameter 
AT C ^ 1 and therefore the line {7 = divides the parame- 
ter space into two parts, the part with dominating CDW 
or SDW instabilities at U > (A" c < 1) and the part 
with dominating superconducting instabilities at U < 
i.e. (A' c > 1). 

At J xy < and U > J xy the spin sector is gapped and 
the spin field gets ordered with vacuum expectation value 
((fi s ) = 7r/2. This leads to the suppression of all instabil- 
ities whose power-law decay is less than 1/r 2 except the 
SDW Z and TS° ones, which show a power-law decay 



(0 S Dw(P)0 SDW .(r)) 
(0 TS o(0)0 TS o(r)} 



r Kc cos(2-Kiyr) 

r -l/Kc 



(51) 
(52) 



Therefore, in the sector A, at J xy < < U the SDW Z is 
the dominating instability in the system and in the sector 
E at J xy < U < the triplet superconducting ordering 
dominates. 

For U < J xy < and for U < J xy /2 the spin sector 
is also gapped but the vacuum expectation value of the 
ordered spin field (ip s ) = 0. This leads to suppression of 
all instabilities with power-law decay less than 1/r 2 except 



Spin gap 
«p s > = ji/2 Kj.<1 

SDW 2 

(TS°) 




Spin gap E 

<(p >=jt/2 K > 1- 












Jxy- 


u 



Luttinger Liquid 

K S >1 K C <1 

. xy 
■ CDW) 



D 

SS 

(C 



«p s >=o 

C CDW (SS + ES) 



Spin gap J *y 

^>=o K>y 

+ KS 
D \V) 



-U=Jxv/2 



Fig. 4. The ground state phase diagram of the 1/4- filled itiner- 
ant t — Jxy — U model. Solid lines indicate borders between the 
weak-coupling limit phases. The thick solid lines denote lines 
in the parameter space where the spin gap opens and indicate 
borders between the weak-coupling limit phases. The thin solid 
line U = marks the transition from the sector of phase dia- 
gram with dominating density instabilities (K c < 1) into the 
sector with dominated superconducting ordering (K c > 1). 
The subleading instabilities are indicated in brackets. 



density-density and singlet superconducting instabilities, 
which show the following power-law decay 

(OcDw(0)O C Dw{r)) r- K - cos(2 ■Kvr) , (53) 
(Oss(O)Oss(r)) ^ (O ES (0)O ES (r)) ~ r" 1 /^ . (54) 

Therefore in the sector C at < U < J xy /2 the CDW or- 
dering dominates, while in the sector D at U < min{0, J xy ] 
the singlet superconducting order is realized. 

Finally, in the sector B ([/ > J xy /2) the LL phase with 
gapless charge and spin excitation spectrum and domi- 
nating easy-plane antifcrromagnctic ordering is realized 
in the ground state. The corresponding correlations show 
a power-law decay 

{O S DW-v(0)OsDwv(r)) - r- K *- l ' K ' , (55) 

where the fixed point value of the spin stiffness parameter 
is given by 

1 



K* = 1 



+ 2^J Jxy ( 2U ~ Jxy) 



(56) 



One can easily show that the CDW and SDW Z cor- 
relations decay faster 

(O C Dw{0)O C Dw(r)) =i 
{O S Dw(0)O S Dw(r)) ^ r~ K °- K ° , 
and are the subleading instabilities in this sector. 



(57) 



3.2.2 v ^ 1/4 

At v 7^ 1/4 the phase diagram is qualitatively the same as 
at quarter-filling (see Fig. 0}. The minor difference con- 
sist of the v dependence of the border lines which mark 
transitions between different phases. In particular: 
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1) The semiplane in the parameter space corresponding 
to dominant density-density correlations (sectors A,B 
and C in the phase diagram) is separated from the 
semiplane with dominant superconducting ordering (sec- 
tors E and D) by the line U = J xy cos(2ttv); 

2) The line U = J xy cos 2 (iris) corresponds to the border 
line between the B and C sectors of the phase diagram. 

To conclude this section, we have shown that the weak- 
coupling phase diagram of the non-half-filled itinerant t — 
J — U model shows a triplet superconducting ordering in 
the sector of the phase diagram with dominating ferro- 
magnetic exchange and easy-plane anisotropy. However, 
as it was shown in Ref. m the case of half-filled band 
and for transverse ferromagnetic exchange stronger then 
some critical value J xy < J xy (—J xy — W, where W is the 
bandwidth), the transition into an insulating phase with 
easy-plane type ferromagnetic ordering takes place. 

Although the transition into the gapless phase with fer- 
romagnetic order (ferrometal) at v ^ 1/2 is a pure finitc- 
bandwith, strong coupling effect the very presence of this 
transition as well as a clear asymmetry between antiferro- 
magnctic (J xy > 0) and ferromagnetic ( J xy < 0) exchange 
can be traced already within the weak-coupling treatment, 
if we take into account the Hartree regularization of the 
hopping amplitude given by £ c ff = t (1 + jJ xy /2Trt). Here 
7(f) > is a band-filling dependent parameter, which 
is of the order of unity for the band-fillings considered 
in this paper. It is clear that in the case of antiferro- 
magnctic transverse exchange the effective bandwidth in- 
creases, while for ferromagnetic coupling it reduces. The 
weak-coupling approach fails when the effective dimen- 
sionlcss coupling constant \g+\ = \J xy \/-KVF becomes of 
the order of unity. Taking into account the Hartree renor- 
malization of the Fermi velocity we estimate the range of 
applicability of the weak-coupling approach as \ J xy \ < irt. 
In this sector of the parameter space with strong easy- 
plane anisotropy \g±\ > \gos\ the soliton mass of the sine- 
Gordon field © is given by M s ~ W exp(-irW/2\J xy \) 
and therefore for \J xy \ < irt the spin gap increases 
with increasing strength of the transverse ferromagnetic 
exchange. However, since in the case of strong ferromag- 
netic exchange J xy < —irt the effective bandwidth W ~ 
t e g tends to zero, the initial increase of the spin gap should 
change into a decrease caused by the collapse of the band- 
width W. As we show below, using the DMRG studies of 
chains up to L = 120 sites, this is indeed the case. For 
Jxy < 0, the spin gap as a function of the parameter J xy 

-irt 

-At. At J < J ( xy 1] the 



shows a bellshaped behavior with maximum at J xy 



and reaches zero at J xy = J xy 



(cl) 



ground state of the system is similar to that of the t — J xy 
model with gapless charge (due to the doping) and gapless 
spin (due to the in-plane, XY character of the exchange) 
excitation spectrum. 



4 Numerical results for J z = and u — 1/4 



we use the density-matrix renormalization-group (DMRG) 
method |37k8l39| . As in our previous study [2] it is ap- 
plied to open chains up to 120 sites keeping typically 400 
states in each block using the infinite-size algorithm to de- 
termine the ground-state properties, including correlation 
functions. 

Below we focus on the case of a quarter-filled band and 
Jz = 0. 
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Fig. 5. Spin gap A s of the itinerant XY model at quarter 
filling, t = 1, J z = and U — (diamonds), [7 = 1 (circles) 
and [7 = 5 (squares). The charge gap A c vanishes for all J xy . 



4.1 Excitation spectrum 

First we determine the low-lying spin and charge excita- 
tions as function of the transverse spin exchange J xy . In 
the numerical calculations, we have found that for open 
boundary conditions the ground state energy for a system 
of L sites with iV-j up-spin and N± down-spin electrons 

Eq L '(N^, N±) remains in the sector with the "z" compo- 
nent of the total spin Sf ot = (i.e. 7V T = TV; = N/2, 
where N = JVj + N± is the total number of electrons) for 
all parameter values studied here. 

As it is commonly used in literature, gaps to exci- 
tations classified as charge excitations are calculated by 
taking the difference between ground-state energies with 
different number of particles. It is convenient to stay in 
the sector with S* t 2 ot = and therefore the charge gap is 
evaluated by 



ML) = \ 



AL) 



A 



N 



N N 
2" _1 '"2 



1 



2Ei L > I '—, '— 



N N 
2~'~2 



.(58) 



In order to check the validity of the picture suggested by 
the bosonization results derived in the previous sections 



We determine the spin gap as the difference between 
the lowest energy in the sector with S^ ot — (7Y| — TYj, = 2) 
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and the ground state energy 

The results for finite chains are extrapolated for L — > 
oo by fitting a polynomial in 1/L. As expected the charge 
gap A c = limi^oo A C (L) vanishes for all values of J xy . 
On the other hand, the spin gap A s = limi,_ >00 A S (L) 
shows a nontrivial behavior where five different regimes 
can be distinguished (Fig. [SJ. In agreement with predic- 




ts -4.0 -3.5 -3.0 -2.5 -2.0 -1.5 -1.0 -0.5 0.0 



Fig. 6. Spin gap A a of the itinerant ferromagnetic XY model 
(J xy < 0) at U = and for fillings v = 1/2, 1/4, 3/8, 1/8. 



tions of the weak-coupling treatment, at U = and weak 
exchange \ J xy \ the spin sector is gapped for both signs 
of the transverse exchange. With increasing antifcrromag- 
netic exchange the spectrum remains gapped, while in the 
case of ferromagnetic exchange, at J xy ~ — 2.5i the spin 
gap starts to decrease and becomes zero at J xy ~ — 3.4i. 
The spin gaplcss phase remains till J xy ~ — 5.8t where the 
spin gap opens once again. 

Thus three gapful regimes (for large ferro- and antifcr- 
romagnetic couplings and for intermediate ferromagnetic 
exchange) are separated by two gapless regimes. As it is 
seen from Fig. a repulsive Coulomb interaction U > 
only changes this behaviour quantitatively. The gap is 
slightly suppressed and the gapful regimes appear at larger 
values of the exchange interaction J xy . 

We have also checked the band-filling dependence of 
the spin-gapped phase for intermediate feromagnetic ex- 
change. As it follows from Fig. |5] the dome-type depen- 
dence of the spin gap on the increasing ferromagnetic 
transverse qualitative remains unchanged at v ^ 1/2. Al- 
though position and value of the maximum of the spin 
gap depends on the band filling (Fig. EJ) the effect of the 
closing of the spin gap at large ferromagnetic exchange is 
band- filling independent. 



0.0 20.0 40.0 60.0 
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HI 

Fig. 7. Correlation functions, plotted against the real space 
distance \i— j\, in the case of strong antiferromagnetic exchange 
Jxy = 3.8t and U = 0. The lower figure shows the decay of the 
correlations, plotted on a double logarithmic scale. Numbers 
in the inset of the lower figure are the exponents. The nota- 
tion "xxx" corresponds to the case of exponentially decaying 
correlations. 



4.2 Correlation functions at U = 

To investigate the nature of ordering in the different phases 
we study the behavior of the correlation functions. In the 
sectors with gapless excitation spectrum we expect the 
usual expression for correlation functions 

C(r) ee (Ot(r)O(O)) ~ A ir ~ 91 + cos (2k F r) A 2 r~ 02 (60) 

consisting of a smooth part decaying with exponent 9\ 
and an oscillating part decaying with 62- In determining 
the asymptotics of correlation functions we focus on the 
dominating part given by 9 = min{#i, #2}- 

In the following we will present results for correlation 
functions in different sectors of the phase diagram. 

4.2.1 Sector I: J xy > 0. 

The SS + ES + CDW phase. 

We start our consideration from the case antiferromag- 
netic exchange. In Fig.[7|we have plotted results of DMRG 
calculations for correlation functions at J xy = 3.8i and 
U = 0. 

As it is clearly seen from Fig. [7| in an excellent agree- 
ment with the bosonization results, the triplet supercon- 
ducting and antiferromagnetic correlations are completely 
suppressed in this case, while the singlet superconducting 
and density-density correlations show an identical power- 
law decay with critical indices equal to one. The very small 
deviation of the numerically evaluated values of the expo- 
nents 9es — 0.97 and Qcdw — 0.98 from the analytically 
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4.2.3 Sector III: Jf y < J xy < Jf y 
The ferrometallic phase 




In Fig. we have shown results of DMRG calculations 
for correlation functions for strong ferromagnetic coupling 
J xy = —5.0t. In this sector, charge and spin gap vanish 
and the system shows properties of a metal with gaplcss 
spin degrees of freedom. As it is clearly seen, the trans- 
verse ferromagnetic and triplet superconducting instabil- 
ities dominate in the ground state. The lower figure in 
Fig. ED displays correlation functions plotted in a double 
logarithmic scale. As it follows from this figure, despite 
the gapless character of spin excitations, the singlet su- 
perconducting correlations are suppressed exponentially. 
At the same time, the SDW Z and CDW correlations also 



Fig. 8. Correlation functions, plotted against the real space 
distance \i — j\ in the case of ferromagnetic exchange J xy — 
—2.8t and U = 0. The lower figure shows the decay of the 
correlations, plotted on a double logarithmic scale. Numbers 
in the inset of the lower figure are the exponents. The nota- 
tion "xxx" corresponds to the case of exponentially decaying 
correlations. 



predicted value 8ss = &es = 8cdw = 1 reflects the high 
accuracy of the obtained numerical results in this sector 
of the phase diagram. 



4.2.2 Sector II: J c x \ < J xy < 0. 



The SDW Z + TS° phase 




Let us now consider the case of ferromagnetic exchange. 
At J xy < and U = 0. the weak-coupling bosonization 
results predict exponential suppression of the CDW and 
singlet correlations, whereas SDW 2 and triplet correla- 
tors TS° show a power- law decay (cf. with Eq. ^ - ^ ). 
Furthermore, they are the dominating instabilities in this 
phase. 

Fig. |H1 displays DMRG results for the correlation func- 
tions. One can clearly observe a strong SDW Z and TS° 
correlation in the ground state. In addition, from the dou- 
ble logarithmic plot one obtains that the in-plane mag- 
netic correlations decay exponentially, while the density 
and singlet-superconducting correlations show an almost 
identical power-law decay at large distances. This is ex- 
pected from the bosonization results, since, due to the 
pinning of the spin Bose field with vacuum expectation 
value (\/2tt K s (p s ) = 7r/2, the in-plane magnetic correla- 
tions which are determined by the dual-field decay expo- 
nentially. The same reason of spin-field pinning leads to 
suppression of the oscillating part in the density correla- 
tions (clearly seen in Fig. |3J) and to the suppression of the 
smooth part of the on-site singlet correlations. The latter 
is the reason why the on-site and extended singlet-pair 
correlations show an identical behavior. 



Fig. 9. Correlation functions, plotted against the real space 
distance in the case of strong antiferromagnetic exchange 

Jxy = —5-0t and U — 0. The lower figure shows decay of the 
correlations, plotted on a double logarithmic scale. Numbers 
in the inset of the lower figure are the exponents. The nota- 
tion "xxx" corresponds to the case of exponentially decaying 
correlations. 



show power-law behavior, but decay slightly faster then 
the in-plane spin correlations. Clearly in this sector the 
system shows very unconventional behaviour, not consis- 
tent with the standard weak-coupling results. 



4.2.4 Sector IV: J xy < Jf y . 

The triplet superconducting phase. 

This is the most intriguing sector of the phase diagram. 
The spin excitation spectrum is gapped. Moreover, as it is 
clearly seen from the Fig. [3] the opening of the spin gap at 
Jxy < J X y — —5.8t is independent of the on-site repulsion 
U. This indicates that in this sector of the phase diagram 
the itinerant nature of the electron system is completely 
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lost and the phase diagram is completely determined by 
the properties of the effective ferromagnetic t— J xy model. 
We expect that this sector of the phase diagram corre- 
sponds to a triplet superconducting phase. 



5 Conclusions 

Motivated by recent experimental findings that show evi- 
dence for the competition or even coexistence of supercon- 
ductivity and magnetism we have continued our studies of 
the ground state properties of an itinerant XY model. Us- 
ing a composite approach based on weak-coupling bosoniza- 
tion and DMRG studies for chains up to L = 120 sites we 
have studied the ground state phase diagram of the itin- 
erant XY model away from half-filling. 



A A, 



Sector II 



Sector III 





Sector I 
SS + ES + CDW 




Fig. 10. Qualitative form of the ground state phase diagram 
for the t — - Jxy — U model at quarter-filling and t = 1 and 
U > 0. The bold red line shows the spin gap as a function of 
the parameter J xy . 



Depending on the model parameters J xy , U > and 
the band-filling v we have found evidence for five different 
phases in the ground state. Fig. 1101 summarizes our find- 
ings for the particular case of a quarter-filled band. Within 
the considered range of parameters the charge gap is al- 
ways zero. The behavior of the spin gap as function of 
the spin-coupling J xy allows to distinguish the following 
different phases: 

1. For antifcrromagnctic interactions J xy > the spin 
gapful metallic phase with dominating singlet super- 
conducting and density-density correlations. 

2. Approaching the line J. xy = the spin gap closes and 
at this point the system shows properties of a Luttinger 
metal. 

3. In the ferromagnetic sector at — 3.8t = J xy < J xy < 
the spin gap is finite. It opens for arbitrary weak fer- 
romagnetic XY-exchange, shows a dome-type shape 
and closes at J xy = J^ y . In this sector of the phase di- 
agram the system the system displays properties of a 
spin gapped metal with coexistence of triplet supercon- 
ducting and spin-density- wave (SDW 2 ) instabilities. 

4. At — 5.8i = J£y < J xy < J£y a ferrometallic phase 
with gapless spin excitations and strongly dominating 
triplet superconducting and transverse ferromagnetic 
instabilities is realized in the ground state. 



5. At J xy < J£ y the spin gap opens once again. We ex- 
pect that in this phase the system shows properties of 
a triplet superconductor. We have to stress that the 
effective ferromagnetic t — J xy model, which governs 
the behaviour of the system in this sector of the phase 
diagram, requires very detailed separate studies. 

Although the presented numerical results arc restricted 
to the case of quarter-filled band, our analysis indicates 
that qualitatively the phase diagram remains similar at 
v 7^ 1/4. However, deviations from the commensurate 
value of the quarter-filled band remove the degeneracy in 
favour of density-density type ordering at J xy cos(27w) < 
and superconducting type ordering at J xy cos(27w) > 0. 

One of the interesting perspective for future studies 
would be the investigation of a genuine ferromagnetic t— J- 
typc model which should help to understand the limit of 
strong Coulomb repulsion better. 
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